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Abstract 

We present a study of orbit jitter and emittance growth 
in a long linac caused by misalignment of quadrupoles. 
First, assuming a FODO lattice, we derive analytical for- 
mulae for the RMS deviation of the orbit and the emittance 
growth caused by random uncorrelated misalignments of 
all quadrupoles. We then consider an alignment algorithm 
based on minimization of BPM readings with a given BPM 
resolution and finite mover steps. 

1 INTRODUCTION 

In this paper we study the emittance dilution of a beam 
caused by quadrupole misalignments in a long linac. To 
suppress the beam break-up instability an energy spread is 
usually introduced in the beam. For the Next Linear Col- 
lider (NLC) [jl[], the rms energy spread within the bunch 
will be of order of 1%. Due to the lattice chromaticity, the 
deflection of the beam by displaced quadrupoles results in 
the dilution of the phase space and the growth of the pro- 
jected emittance. 

The effect of lattice misalignments has been previously 
studied in many papers. A qualitative analysis and main 
scalings were obtained in Ref. and detailed studies 
with intensive computer simulations are described in Refs. 

^ H]. The purpose of this paper is to develop a simple 
model based on a FODO lattice approximation for the linac 
which allows an analytic calculation of the emittance dilu- 
tion. The model can be also generalized, to include a slow 
variation of the lattice parameters, as well as variation of 
both beam energy and the energy spread [^. 

Throughout this paper we assume that the number of 
quadrupoles in the linac is large, N ^ 1, and neglect terms 
of the relative order of N^^ in the calculations. For future 
linear colliders with the center of mass energy in the range 
of 1 TeV, typically N ~ 10'^, and N^^ is indeed a small 
number 

2 BEAM ORBIT IN MISALIGNED 
LATTICE 



Let us consider a FODO lattice with a cell length / and a 
phase advance fj, per cell, consisting of N thin quadrupoles 
as shown in Fig. ^ The focal length of the quadrupoles is 
equal to ±F where the positive and negative values of F 
refer to the focusing and defocusing quadrupoles respec- 
tively. The beam is injected in the linac at the center of the 
first quadrupole, at s = 0, with the zero offset and angle, 
and the beam emittance is measured at the center of the 
last, A^th, quadrupole. For the beam position (horizontal or 
vertical) at the locations of the quadrupoles we will use the 
notation xi,X2, ■ ■ ■ ,xn-i,xn, and the orbit angle at the 
center of the fcth quadrupole is denoted by x'f. . The initial 
conditions for the orbit are xi = x'l = 0. Note that due to 
our choice of positions, the derivative of the beta function, 
and hence the Twiss parameter a, at all locations 1 through 
N, are equal to zero. 

We now assume that each quadrupole in the lattice is 
misaligned in the transverse direction relative to the axis 
of the linac by ^i, (1 < i < N), where are random, 
uncorrelated numbers. Due to the deflection by misaligned 
quadrupoles, the original straight orbit will be perturbed. 
The offset Xi can be found as 

Xi = ^RikOk, (1) 

k<i 

where Ri^ is the (1,2) element of the transfer matrix R 
and 6i is the deflection angle resulting from the offset of 
the ith quadrupole, di = ±f ^ / F, for the focusing and de- 
focusing quadrupoles. We have = y/ Pi/3k sin Ai/^ik, 
where the betatron phase advance Aipik between fcth and 
ith quadrupoles (k < i) is A-0ife = (1/2) (« — fc)/i. 

We will also need the orbit angles x'^ where the prime 
denotes the derivative with respect to the longitudinal co- 
ordinate s. For x[ we have 

k<i 

where Git, is the (2, 2) element of the transfer matrix, 
Gik = \/ hi Pi cos A^pik, (note that, due to our choice, 
a, = 0). 
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Figure 1 : FODO lattice of a linac. Beam positions are mea- 
sured at the center of each quadrupole. 
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3 RMS VALUE FOR THE BEAM OFFSET 

To characterize the deviation of the orbit from the linac 
axis, we will calculate the average value {x%), where the 
angular brackets denote averaging over all possible values 
of ^. We assume that the average offset {^i) vanishes, hence 

(xn) = 0. 

For the lattice shown in Fig. |l] the deflection angle Ok 
due to the misaligned fcth quadrupole is given by 9k = 



{—Ip^k/F, and the beam offset at the end of the linac is 



XN 



(3) 



A;=l 



For the variance of xn we have 



F2 



Af-1 

1=1 



2 

Nh 



(4) 



where we have used {^k£,i) = (C^)'5fci, with (f^) being the 
variance of the random variables ^i. To calculate the sum 
in Eq. (Q), one can average Rj^^ over the betatron phase 
\PnPi- One finds, 



value 



'N 



4F2 



-(/3max + Anin) - tan |. (5) 



We see that the rms value (xj^)^^"^ scales as N^^"^, which 
is a characteristic feature of the random walk motion. 

In a similar fashion, one can find the rms angular spread 
orbits {x'j^)^^^ at the end of the linac. Starting from the 
general expression 



N 



fe=l 



(6) 



and performing the same averaging as for derivation of Eq. 
(||), one finds {x'%) — {xj^)/Pj^, where {xj^) is given by 
Eq. (|). 



4 CHROMATIC EMITTANCE GROWTH 

When the beam has a nonzero energy spread, due to the 
chromaticity of the lattice, the misalignments cause an ef- 
fective emittance growth of the beam [||. We will calculate 
the emittance increase, assuming that the beam energy E 
and the relative energy spread in the beam S are constant. 
We will also assume that the resulting emittance growth is 
much smaller than the initial emittance of the beam. In this 
case, we can use the following formula for the final emit- 
tance growth 



Ae 



1 



+ /3iv((Ax' 



{Ax')f 



(7) 



where Ax and Ax' are the spread in the coordinate and 
the angle within the bunch at the and of the linac, and the 
angular brackets with the subscript denote a double av- 
eraging: first, averaging over the random misalignment of 
the quadrupoles and then averaging over the energy dis- 
tribution function in the beam. We will assume that the 
energy spread in the beam S is so small, that one can use a 
linear approximation for calculation of Aa; and Ax', Ax = 
6 ■ xns = S ■ dxN/dS and Ax' = S ■ x'^^ = 6 ■ dx'^/dS. 



Since (^j) = 0, hence {Ax) = {Ax') — 0. In this approxi- 
mation Eq. reduces to 



N \ ^ NS 



(8) 



where 6'^ is the variance of the energy spread within the 
beam. 

To calculate X]^s and x'j^^g we need to take the derivatives 
of Eqs. (|]) and ^ with respect to S. For a long linac, the 
dominant contribution to Ae comes from the dependence 
of the phase advance Atpik versus energy, so we need to 
differentiate only sin Aipn- (or cos Atpik) terms in the sum. 
Calculation gives 

P],'{x%s)-pN{x'ls)-lN'^-^tan'^. (9) 
which gives for the emittance dilution 



Ae= W^tan^^, 
3 I 2 



(10) 



As we see, the increase in the emittance scales with the 
number of quadrupoles as N^. 

In the above derivation, to find the dispersion of the beam 
at the end of the linac, we explicitly differentiated Eq. ^ 
with respect to the energy. One can use another formula for 
computing 9xAr/9(5 [0], 



dx 



N-l 



N 
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(11) 



k=l 



that takes into account that the dispersion is generated due 
to the offset of the particle relative to the center of the 
quadrupole, and propagates downstream with the same ma- 
trix element RMk- 

5 VERY LONG LINAC 

Increasing the length of the linac and the number of 
quadrupoles brings us to the regime where Eq. (|lO|) is 
not valid any more. The transition occurs when the phase 
advance over the length of the linac due to the energy vari- 
ation (5 becomes comparable to 7r/2, NS-dfi/d6 ^ tt/2. In 
this case, the differential approximation Aa; = S ■ dxN /d5 
that was used in Section Q is not valid any more, and the 
scaling Ae cx breaks down. 

We can estimate the emittance dilution in this regime, us- 
ing the following arguments. Let us denote by Ic the deco- 
herence length in the linac such that {lc/l)5 ■ dji/dS ^ -k 12 
(I is the FODO cell length). When the beam passes the dis- 
tance Ic, due to filamentation, the betatron oscillations of 
the beam are converted into the increased emittance, and 
the subsequent motion becomes uncorrelated with the pre- 
viously excited betatron oscillations. The emittance growth 
on the distance 1^ is given by Eq. ( 10), in which iV — 2lc/l, 



(a 
I 



tan' 



3 M 



{i 



^- (12) 



The total emittance increase in the Hnac of length I l in this 
regime is equal to Aec multiplied by the number of coher- 
ent distances Il/Ic in the linac 

Ae = A.>^M!L^M!l. (13) 
Ic P tan ^ I tan ^ 

Note that if the hnac length < Ic, the emittance dilu- 
tion is reversible in principle - the initial beam emittance 
can be recovered by taking out the dispersion generated by 
the misaligned quadrupoles downstream of the linac. For 
very long linacs, when 1^ > Ic, the emittance growth be- 
comes irreversible due to the phase space filamentation. 

6 ALIGNMENT WITH ACCOUNT OF 
BPM ERRORS AND FINITE MOVER 
STEPS 

Measuring the beam position at each quadrupole, with the 
knowledge of the lattice functions, allows us to find the 
quadrupole offsets ^i. Moving the quadrupoles by distance 
— would position them in the original state, and restore 
the ideal lattice. Of course, in reality, there are many fac- 
tors, such as wakefields and measurement errors, that do 
not allow to perfectly align the lattice. Here we will study 
two such effects: errors associated with the BPM measure- 
ments, and finite step of the quadrupole movers. 

Consider first the effect of BPM errors. Due to the finite 
resolution of BPMs the measured vector of the beam trans- 
verse offsets = {xf^ . . . xfj ) differs from the exact 
values X = {xi . . . xn) by an error vector e, X*^ = X+e, 
where e = (ei . . . ejv)- The errors are small relative to the 
measured values, \ei\ ^ \xi\. We assume that the BPMs 
are built in the quadrupoles, and the quadrupole displace- 
ment ^fe also moves the center line of the BMP, so that BPM 
reading is x^^ = Xk — £,k + Sfc. Using the measured offsets 
x^^ we infer the quadrupole offsets (k from the following 
equation 

+ C.-Ei?.fe(-lf|. (14) 
fe=i 

Note that without errors, ek — 0, we would find from Eq. 
( p4| ) the correct value (^k = Cfc- Measurement errors ek 
cause the inferred values of the offsets differ from the true 
ones, Ck^^k- 

We then ahgn the lattice by moving the quadrupoles by 
distance —(k- After the alignment the corrected beam orbit 
Xi does not vanish: 

fe=i 

= X, -xf -G = -e.+e»-0- (15) 

Since the quadrupoles after alignment are located at ^fc— Cfc, 
the beam offset relative to the center of the quadrupole, 
Xk — (ffe — Cfc). is equal to — e^. This allows us to use 



Eq. ( pi] ) to find the emittance dilution in the linac after the 
alignment, 

N-l 

xjv5--5I^JVfe(-l)''"'y. (16) 

A:=0 

Assuming that are uncorrected random numbers makes 
the problem equivalent to the orbit equation with the 
result given by Eq. (Q), 

/3^'(2;L-)=47V^tan|. (17) 

We see that the rms value of the dispersion at the and of 
the linac after alignment scales as \/~N. Calculating in a 
similar way the variance of the derivative x^^, gives the 
chromatic emittance growth after alignment, 

Ae 4iVF-^ tan I . (18) 

Let us now assume that in addition to the BPM errors the 
quadrupole movers have a finite step so that the final posi- 
tion of the quadrupoles Qk after alignment is ^k — Cfc + '''k, 
where as above, C,k is the offset inferred from the mea- 
surements (and containing BPM errors), and Vk is the 
quadrupole movement error. Again, we assume that are 
random, uncorrected numbers, and of course uncorrected 
with the BPM errors e^. For the beam orbit after alignment 
we now have 

i-l 

Xi = -ei + - Ci + ^ RikTk (19) 
fe=i 

with the resulting emittance growth that is a combination 
ofEqs. and (0), 

Ae = 47V^^ tan ^ + -^iV^ tan^ ^ . (20) 

From this equation, it follows that for a large N, the con- 
tribution of the movers errors becomes more important and 
imposes tighter tolerances on the movers. 
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